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ABSTRACT 

The  sudden  onset  of  the  northern  summer  monsoon  has  led  to 
many  theories,  as  to  the  cause,  one  of  these  being  a  possible 
relation  to  the  elevated  Tibetan  highland  that  dominates  the 
area. 

Using  a  zonally  symmetric  numerical  model  driven  by  a 
specified  equilibrium  heating  function,  two  experiments  are 
conducted.   The  first  experiment  integrates  the  model  over  a 
180-day  period  on  a  flat  earth,  linearly  changing  the  equilib- 
rium temperature  from  January  to  July.   The  second  experiment 
repeats  the  first  except  for  the  inclusion  of  the  topography 
of  85°E  longitude. 

In  both  experiments,  actual  velocities  are  underestimated 
with  the  topography  experiment  providing  the  better  results. 
Development  of  the  monsoon  was  poor  in  both  experiments,  due 
mainly  to  the  inability  of  the  model's  vertical  resolution  to 
simulate  the  extreme  change  in  topography  elevation  and  the 
unrealistic  approach  of  a  linear  temperature  change.   Possible 
improvements  for  future  studies  are  suggested. 
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I,   INTRODUCTION 

Monsoons,  meaning  "season"  from  the  Arabic,  are  complex, 
seasonally  dependent  weather  systems  covering  the  entire 
tropics  and  sub  tropics.   The  northern  summer  monsoon  begins 
in  June  and  is  characterized  by  a  surface  monsoon  trough 
over  India,  with  an  upper  level  anticyclone  over  the  Tibetan 
plateau.   To  the  south  of  this  Tibetan  anticyclone,  near  10°N 
to  15°N,  is  a  strong  easterly  jet  stream  which  encompasses 
almost  the  entire  tropics,  the  maximum  of  the  jet  stream 
velocities  being  centered  at  approximately  85°E  longitude. 

Of  interest  to  meteorologists  for  many  years  has  been  the 
rapidity  with  which  the  onset  of  the  monsoon  activity  occurs. 
During  the  winter  months,  the  Asia  continent  is  under  the 
influence  of  a  subtropical  westerly  jet  stream  at  the  upper 
levels.   Over  a  very  brief  period  of  time,  measured  in  days, 
the  westerly  jet  stream  migrates  toward  more  northern  lati- 
tudes and  is  replaced  in  the  tropics  by  the  tropical  easterly 
jet  stream.   Coupled  with  this  replacement  is  the  surface 
monsoon  trough  development  and  its  associated  surface  south- 
westerly winds,  winds  laden  with  moisture  from  the  oceans. 
These  southwesterly  winds  cause  widespread  rainfall  as  they 
reach  the  southern  Asia  continent. 

The  need  to  understand  the  sudden  onset  of  the  summer  mon- 
soon has  led  to  many  theories  among  meteorologists.   Yin  (1949) 
has  proposed  that  this  onset,  as  represented  by  the  sudden 
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replacement  of  the  westerly  jet  by  the.  easterly  jet  in  the 
tropics,  may  be  related  in  some  way  to  the  effect  of  the 
Tibetan  highlands,  an  elevation  that  is  very  striking  in 
nature. 

The  purpose  of  this  study  is  to  examine  such  a  possibility 
by  using  a  numerical  model  to  simulate  the  development  of  the 
summer  monsoon  with  different  experiments.   For  the  first 
experiment,  a  zonally  symmetric  January  steady  state  condition 
on  a  flat  earth  is  obtained.   Starting  from  this  condition, 
the  model  is  integrated  for  180  days  with  the  mean  thermal 
condition  changing  from  January  to  July.   The  final  state 
represents  a  July  steady  state  condition.   In  the  second  ex- 
periment, the  distribution  of  topography  along  85°E,  which  is 
near  the  center  longitude  of  the  Tibetan  highland,  is  included 
in  the  model.   The  180-day  integration  is  repeated. 

It  is  hoped  that  by  comparing  the  results  of  these  two 
experiments,  some  understanding  of  topographical  effects  upon 
the  evolution  of  events  that  leads  to  the  development  of  the 
summer  monsoon  may  be  obtained. 


II,   MODEL  AND  PROCEDURE 

A.   BASIC  MODEL 

In  this  study,  the  model  formulated  by  Monaco  and  Williams 
[1975]  is  used.   This  model  uses  the  primitive  equations  in 
spherical  coordinates  with  "sigma"  as  the  vertical  coordinate. 
The  horizontal  distribution  of  variables  is  on  a  4-degree 
resolution.   The  model  is  made  zonally  symmetric  by  removing 
all  zonal  variations.   The  troposphere  is  divided  into  seven 
equally  spaced  "sigma"  levels  with  three  main  reporting  levels 
The  sigma  coordinate  is  defined  as: 

P  "  Pt 
o    =   =-   ,  (2.1) 

where  p  is  the  pressure  of  the  sigma  level,  p   is  the  height 
of  the  constant  tropopause  which  is  set  at  100  mb ,  and  tt  is 
the  terrain  pressure  defined  as: 

*  =  Ps  -  Pt  ,  (2.2) 

where  p   is  the  surface  pressure. 

Time  integration  is  comprised  of  continuous  sections  of 
one  Matsuno  and  four  leapfrog  steps  with  time  increments  of 
6  minutes.   All  data  for  this  study  are  interpolated  to  con- 
stant pressure  surfaces  of  250  mb,  500  mb ,  and  850  mb . 
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B.   SPECIFICATION  OF  HEATING  AND  TOPOGRAPHY 

The  model  is.  driven  by  a  heating  function,  H,  which  con- 
tinuously adjusts  the  model  atmosphere  temperature  to  a 
specified  equilibrium  temperature,  T* ,  using  the  following 
formula: 

H  =  -  [  T  ~  T*Cy,<r)  j  f  C2.33 

where  T  is  the  model  atmosphere  temperature  and  t  is  a  speci- 
fied adjustment  time,  set  at  two  days  for  this  study. 

In  the  no- topography  experiment,  the  equilibrium  tempera- 
tures for  January  and  July  are  extracted  from  observations  of 
Palmen  and  Newton  [1969]  for  the  January  and  July  climatology 
(Figures  1  and  2,  respectively).   The  dashed  lines  indicate 
the  horizontal  sigma  levels  used  in  obtaining  the  tempera- 
tures, in  this  case  assuming  the  model  levels  to  be  250  mb , 
550  mb,  and  850  mb . 

Because  the  temperatures  in  July  at  85°E  are  considerably 
warmer  than  the  zonally  averaged  temperatures,  the  equilibrium 
temperature  for  July  is  modified  by  the  addition  of  a  departure 
from  the  zonal  mean,  at  85°E,  as  observed  by  Krishnamurti  [1971] 
and  shown  in  Figure  3.   At  the  250  mb  level,  the  full  depar- 
ture is  added  to  the  mean  temperature.   At  the  550  mb  level, 
50  percent  of  the  departure  is  added.   At  the  850  mb  level, 
20  percent  of  the  departure  is  added  to  the  mean  temperature. 

In  the  topography  experiment,  topography  as  a  function  of 

latitude  is  shown  in  Figure  4.   It  is  introduced  through  a 

gradual  change  in  the  surface  geopotential ,  cj>  ,  as  given  by 

the  following  formula: 
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Figure  3:   Deviations  from  the  zonal  mean  July 
temperature  at  85°E  longitude. 
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Figure  4:   Topography  elevation  at  85°E  longitude 
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<Ds  =  4)f  (  sin2  (  t^))  •  (2.4) 

Here,  4>.p  is  the  specified  final  surface  geopotential  value; 
t   is  the  total  amount  of  time  required  for  mountain  building, 
specified  as  24  hrs  for  this  study;  and  t   is  the  number  of 
incremental  changes  to  occur  within  the  mountain  building 
item,  specified  as  48  hrs  for  this  experiment. 

The  gradual  build-up  of  the  topography  is  to  avoid  the 
undesirable  numerical  shock  of  a  sudden  change  invariable 
magnitudes.   The  same  January  equilibrium  temperature  for  the 
no- topography  experiment  is  used  in  this  experiment  before 
the  introduction  of  topography,  but  after  the  topography  is 
fully  developed  the  temperature  specified  along  the  dotted 
lines  in  Figure  1  is  used.   These  curves,  which  are  initial 
sigma  levels,  are  patterned  after  the  topography. 

The  July  equilibrium  temperature  for  this  experiment  is 
extracted  from  Figure  2,  following  the  same  dotted  lines. 
Modif ication  of  the  July  temperature  is  again  accomplished 
by  adding  the  same  deviation  from  the  mean  as  used  in  the 
no- topography  experiment. 

C.   NO- TOPOGRAPHY  EXPERIMENT 

In  this  experiment,  pressure  levels  are  initialized  at 
standard  atmospheric  values  with  all  velocities  being  zero. 
The  initial  temperature  is  specified  as  the  January  equilibrium 
temperature  as  shown  in  Figure  1.   The  model  is  started  and 
allowed  to  run  for  a  period  of  90  days,  sufficiently  long  to 
establish  a  quasi-steady  state  condition. 
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The  integration  is  continued  from  this,  steady  state  con- 
dition and  allowed  to  run  for  a  period  of  another  180  days. 
Throughout  this  run,  the  equilibrium  temperature  is  continu- 
ously changed,  by  a  simple  linear  interpolation,  to  reach  the 
equilibrium  temperature  of  July. 

At  the  end  of  the  180-day  run,  the  integration  is  allowed 
to  continue  for  an  additional  90  days  with  the  July  equilibrium 
temperature  remaining  fixed  in  order  to  achieve  a  quasi-steady 
state  condition  for  July. 

There  is  no  significant  difference  between  the  circulation 
patterns  at  the  end  of  the  180-day  run  and  those  at  the  end  of 
the  additional  90  days  using  constant  July  temperatures.   Thus, 
it  is  apparent  that  the  model  adjusts  well  to  the  equilibrium 
heating  function  and  any  day  within  a  particular  simulation 
period  may  be  considered  to  represent  an  instantaneous  quasi- 
steady  state  condition. 

D.   TOPOGRAPHY  EXPERIMENT 

The  procedure  for  this  part  of  the  experiment  is  similar 
to  the  no- topography  case  except  for  the  introduction  of  to- 
pography.  To  compensate  for  the  mass  changes  due  to  the  volume 
of  atmosphere  change,  the  surface  pressure  is  specified  as 
860  mb ,  a  figure  experimentally  determined  to  produce  an  over 
ocean  surface  pressure  of  1015  mb . 

The  same  January  equilibrium  temperature  used  for  the  no- 
topography  experiment  is  held  fixed  and  the  model  is  integrated 
for  a  period  of  45  days,  allowing  the  model  to  adjust  to  the 
artificially  low  initial  surface  pressure.   At  this  point, 
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the  topography  is  introduced  as  previously  discussed  and  the 
equilibrium  temperature  is  readjusted  to  conform  with  the 
topography  sigma  levels  as  shown  by  the  dotted  line  in  Figure 
1.   An  additional  45  days  is  allowed  after  the  introduction 
of  topography  and  January  temperature  change  to  achieve  a 
January  steady  state. 

From  this  steady  condition,  the  integration  is  allowed  to 
proceed  for  180  days,  again  allowing  the  equilibrium  tempera- 
ture to  change  linearly  from  January  to  July,  the  July  temper 
ature  being  obtained  as  previously  discussed. 

The  July  equilibrium  temperature  is  again  held  fixed  and 
the  model  is  allowed  to  continue  for  an  additional  90  days. 
There  is  again  no  significant  difference  between  the  wind 
patterns  at  the  end  of  the  180-day  integration  and  the 
additional  90-day  run. 
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III.   RESULTS 

The  principal  circulation  features  of  interest  in  this 
study  are  the  Tibetan  high  and  the  tropical  easterly  jet 
stream.   The  simulation  of  the  development  of  these  features 
by  the  model  are  compared  with  observations  at  85°E,  mainly 
those  of  Sadler  [1975] . 

A.   SIMULATED  CIRCULATIONS  IN  JANUARY  AND  JULY 
1 .   No-Topography  Experiment 

The  January  steady  state  condition  without  the  inclu- 
sion of  topography  is  represented  by  the  zonal  wind  as  a 
function  of  latitude  and  elevation  in  Figure  5.   The  zonal 
wind  at  250  mb  as  a  function  of  latitude  is  compared  with 
that  observed  by  Sadler  at  200  mb  in  Figure  6.   The  observed 
westerly  flow  occupies  the  region  north  of  approximately  14°N 
with  a  well  defined  jet  core  at  32°N.   The  maximum  velocity 
in  this  core  is,  on  the  average,  44  m  sec   .   The  model  pro- 
duces a  broad  and  weak  westerly  flow  to  northward  of  4°N. 
There  is  a  definite  lack  of  a  well  defined,  sharp  jet  maximum, 
the  maximum  velocity  being  21.6  m  sec    at  40°N.   At  all 
latitudes,  the  velocities  are  grossly  underdeveloped. 

The  July  steady  state  condition  is  represented  by  the 
zonal  wind  as  a  function  of  latitude  and  elevation  in  Figure 
7.   The  zonal  wind  at  250  mb  as  a  function  of  latitude  is 
compared  with  Sadler's  observation  at  200  mb  in  Figure  8. 
The  replacement  of  the  westerly  jet  stream  by  the  tropical 
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Figure  6:   Simulated  zonal  wind  at  250  mb  vs.  observed  zonal 
wind  at  200  mb  for  January,  no  topography. 
Units  are  m  sec~l. 
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Figure  8:   Simulated  zonal  wind  at  250  mb  vs.  observed 

zonal  wind  at  200  mb  for  July,  no  topography. 
Units  are  m  sec"-'-. 
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easterly  jet  is  quite  evident  in  both  data.   The  observed 
easterly  jet  axis  is  centered  at  12°N  with  a  maximum  velocity 
of  22  m  sec"  .   In  this  experiment,  the  model  produced  a  more 
realistic  result  than  the  January  state.   The  simulated 
tropical  easterly  jet  is  clearly  developed  with  a  core  axis 
at  16°N.   The  maximum  velocity  in  the  core  is  23  m  sec   ,  in 
close  agreement  to  the  observed  value. 
2 .   Topography  Experiment 

The  January  steady  state  condition  with  the  inclusion 
of  the  topography  is  displayed  in  Figure  9  which  gives  the 
zonal  wind  as  a  function  of  latitude  and  elevation.   The  zonal 
wind  at  250  mb  is  again  compared  with  Sadler's  observation  at 
200  mb  in  Figure  10.   The  underestimate  of  winds,  a  problem 
in  the  no- topography  case,  is  still  evident  but  not  to  such 
an  extreme  degree.   For  example,  at  32°N,  the  no- topography 
value  of  16  m  sec    is  increased  to  26  m  sec   .   It  is  still 
below  the  44  m  sec   observation,  however.   The  simulated 
pattern  changes  markedly  from  the  no- topography  case.   The 
simulation  shows  a  clearly  bimodal  distribution  with  jet  stream 
cores  of  25  m  sec    at  20°N  and  26  m  sec    at  32°N,  the  latter 
coincides  with  the  observed  latitude  of  the  westerly  jet  core. 
The  latitude  band  between  these  two  jet  cores  contains  the 
vertical  face  of  the  model  topography. 

The  July  steady  state  condition  is  shown  in  Figure  11 
by  the  zonal  wind  as  a  function  of  latitude  and  elevation. 
Figure  12  gives  the  comparison  of  the  simulated  zonal  wind 
at  250  mb  and  the  observed  zonal  wind  at  200  mb  bv  Sadler  as 
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Figure  10: 


Simulated  zonal  wind  at  250  mb  vs.  observed 
zonal  wind  at  200  mb  for  January  with  topography 


Units  are  m  sec 
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Figure  12:   Simulated  zonal  wind  at  250  mb  vs.  observed 

zonal  wind  at  200  mb  for  July,  with  topography 
Units  are  m  sec~l. 
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a  function  of  latitude.   In  this  case,  the  near  match  in  the 
velocities  in  the  latitudes  of  the  tropical  easterly  jet  has 
diminished  somewhat.   The  jet  core  is  shifted  north  to  20°N 
and  the  velocity  is  27  m  sec   ,  an  overestimate  from  obser- 
vation by  approximately  7  m  sec 

There  are  no  good  observational  data  at  85°E  longitude 
with  which  to  compare  the  500  mb  and  850  mb  flow  simulations.' 
The  patterns  are  consistent  with  the  250  mb  patterns,  however, 
as  to  placement  of  features  and  intensities  of  circulations. 

B.   DEVELOPMENT  OF  THE  MONSOON  FLOW 

The  zonal  wind  at  250  mb  produced  in  the  no- topography 
experiment  is  shown  as  a  function  of  latitude  and  time  in 
Figure  15.   The  zonal  wind  at  250  mb  produced  in  the  topography 
experiment  is  shown  as  a  function  of  latitude  and  time  in 
Figure  14. 

In  both  the  topography  and  no- topography  experiments,  the 
seasonal  reversal  in  the  zonal  wind  at  the  250  mb  level  is 
completed  gradually  over  the  six-month  period. 

The  tropical  easterly  jet  first  appears  at  approximately 
12°N  in  mid-March  for  the  no- topography  experiment.   In  the 
topography  experiment,  it  is  developed  at  approximately  20°N 
in  early  March,  thus  a  shift  of  the  jet  axis  to  the  north  by 
8  degrees  of  latitude  with  the  introduction  of  topography  is 
produced.   This  displacement  in  the  core  is  not  observed  and 
thus  must  be  considered  incorrect.   In  the  real  atmosphere, 
the  appearance  of  the  tropical  easterly  jet  actually  begins 
in  mid-May  and  is  centered  at  approximately  12°N. 
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Because  of  the  early  arrival  of  the  easterly  jet  in  both 
experiments,  the  simulated  timing  of  the  monsoon  development 
is  very  poor,  with  the  no- topography  case  being,  paradoxically, 
the  better  of  the  two,  however. 

The  gradual  acceleration  of  the  zonal  winds  for  both  cases 
is  apparently  related  to  the  specified  linear  change  of  the 
equilibrium  temperatures.   In  the  real  atmosphere,  the  temper- 
ature change,  i.e.,  the  reversal  of  north-south  temperature 
gradient,  occurs  rather  suddenly  in  mid-May,  so  that  the 
simulation  of  the  change  in  temperature  is  unrealistic  and 
it  causes  the  gradual  reversal  of  the  zonal  flow. 
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IV,   DISCUSSION 

The  results  of  the  two  experiments  indicate  that  the 
simplifications  used  in  the  numerical  model  are  unrealistic. 
These  simplifications  can  be  broken  into  two  general  areas, 
one  related  to  the  model  itself  and  a  second  related  to  the 
data  used. 

A.  MODEL  SIMPLIFICATIONS 

For  both  experiments,  a  zonal  symmetry  assumption  was 
imposed.   It  appears  that  an  exclusion  of  zonally  asymmetric 
eddy  motions  may  lead  to  an  unacceptable  circulation  pattern, 
due  to  the  importance  of  nonlinear  processes,  as  shown  by 
Krishnamurti  [1971]  and  others.   A  more  severe  problem  of  the 
model  is  perhaps  the  vertical  resolution,  as  evident  by  the 
extremely  poor  simulation  of  the  time  evolution  of  the  250  mb 
zonal  flow  in  the  topography  experiment.   Only  three  levels 
in  the  troposphere  are  apparently  insufficient  because  of  the 
steepness  of  the  Himalaya  mountains.   As  a  result,  the  inter- 
polation of  the  simulated  data  from  sigma  levels  to  constant 
pressure  surfaces  suffers  from  the  crude  vertical  resolution. 

B.  DATA  SIMPLIFICATIONS 

Clearly  shown  in  the  model  results,  the  zonal  velocity  is 
underestimated  in  January  for  both  the  topography  and  no- 
topography  experiment,  although  the  July  velocities  are  in 
somewhat  better  agreement  with  observation.   As  stated 
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previously,  the  July  equilibrium  temperature  specified  in  the 
model  has  been  modified  from  the  mean  by  observations  of 
Krishnamurti  [1971] .   No  similar  modification  was  made  in 
January  as  only  the  zonal  mean  temperatures  are  used.   Because 
the  wind  is  quite  geostrophic  and  is  coupled  with  temperature 
gradient  through  the  thermal  wind  equation,  an  error  in  tem- 
perature will  precipitate  an  error  in  winds.   It  is  therefore 
necessary  to  improve  the  specification  of  the  thermal  field 
for  January. 

The  specification  of  seasonal  temperature  changes  by  the 
linear  interpolation  between  January  and  July  produces  the 
obvious  and  unrealistic  gradual  and  almost  linear  wind  changes. 
It  appears  necessary  that  the  time  series  of  the  actual  temper- 
ature data  must  be  used  to  adjust  the  thermal  state  of  the 
model  atmosphere. 

Only  when  this  realistic  temperature  change  and  an  adequate 
vertical  resolution  are  used  in  the  model  can  the  mechanical 
effect  of  the  mountains  on  the  zonally  symmetric  circulation 
be  adequately  studied. 
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